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1.0  INTRODUCTION 


The  thought  of  possible  flights  of  supersonic  transports  (SST's)  has 
made  it  increasingly  necessary  to  ascertain  accurately  and  reliably  the  lateral 
extent  of  sonic  boom  on  either  side  of  the  flight  track  -  not  necessarily  the 
same  -  in  the  presence  of  wind  and  temperature  gradients.  This  type  of  infor¬ 
mation  is  necessary  in  planning  the  location  of  possible  SST  routes  so  that 
heavily  populated  areas  are  outside  the  sonic  boom  corridor.  The  width  of  the 
sonic  boom  corridor  is  limited  by  the  normal  atmospheric  temperature  gradients 
which  bend  or  refract  the  sound  rays;  outside  a  certain  width  these  rays  do  not 
reach  the  ground,  thus  producing  a  sound  ’shadow’  from  there  on  out.  Wind 
patterns  modify  the  refraction  of  the  rays.  Calculation  of  the  corridor  width 
(and  lateral  shift  in  case  of  a  side  wind)  depend  heavily  on  the  theory  of 
geometrical  acoustics. 

Classical  theories  of  geometrical  acoustics  have  been  developed  in 
detail  in  Refs.  1,  and  2.  A  detailed  analysis,  and  to  some  extent  synthesis 
of  established  theories  for  sonic  boom  propagation  in  a  horizontally  stratified 
atmosphere  with  winds  was  presented  in  Bef.  3*  The  analysis,  though  similar  in 
some  respects  to  the  present  calculations,  is  far  too  complicated  for  easy  appli¬ 
cation  and  the  emphasis  in  that  report  is  on  pressure  signatures  rather  than 
corridors . 


Kane,  et  al.  Ref.  5,  applied  the  theory  of  Ref.  k  to  predict  the 
variations  in  overpressure  which  would  occur  on  the  ground  as  a  result  of  the 
variations  in  the  atmospheric  properties  between  the  airplane  and  the  ground. 
Although  considerable  information  is  available  in  this  report,  it  is  difficult 
to  extract  it  because  in  many  cases,  the  parameters  for  which  the  curves  were 
plotted  were  missing. 

Calculations  similar  to  those  of  Ref.  5  but  without  winds  were  made 
earlier  by  Randall,  Ref.  6.  Dressier,  et  al,  Ref.  7  calculated  ray-ground 
intersections  in  the  presence  of  wind  and  temperature  gradients,  using  Snell's 
law.  The  method  involves  slicing  the  medium  into  piecewise-constant  thin  layers 
for  temperature,  wind  magnitude  and  direction,  and  applying  Snell’s  law  for  the 
refractive  behaviour  at  the  discontinuities  at  each  interface.  Thus,  the  method 
employs  straight  line  rays  and  constant  wind  magnitude  and  direction  witnin  each 
layer.  While  the  technique  is  basically  correct,  it  is  less  realistic  and  less 
accurate  than  a  solution  employing  linear  variation  of  vector  wind  and  tempera¬ 
tures  between  given  points.  Reed,  et  al.  Ref.  8,  used  an  acoustic  ray  tracing 
procedure  to  determine  the  ground-level  sound  patterns  for  straight,  level 
supersonic  flight  in  arbitrary  weather  conditions.  The  authors  however  treated 
the  ray  as  the  orthogonal  trajectory  of  the  wave  front  -  a  procedure  which  would 
be  correct  in  the  absence  of  winds.  For  small  wind  velocities,  the  error  involved 
in  this  procedure  is  small,  but  for  moderate  and  strong  winds  the  procedure  leads 
to  totally  unacceptable  results. 

It  is  thus  desirable  that  a  single,  accurate,  set  of  calculations  of 
sonic  boom  corridors  in  the  presence  of  wind  and  temperature  gradients  for 
aircraft  in  steady  level  flight  be  carried  out  with  the  aim  of  providing  prac¬ 
tical  information  needed  for  planning  the  possible  location  of  SST  routes.  Such 
calculation  is  the  prime  objective  of  this  study.  In  the  past,  most  investigators 
have  used  the  conmonents  of  wind  velocity  along  the  aircraft  axes  as  the  variable 
wind  parameters.  Since  however,  the  instrument  panel  of  the  aircraft  will 
presumably  provide  such  data  a s  aircraft  heading,  Mach  number,  altitude,  wind 
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magnitude  and  direction,  it  serves  a  useful  practical  purpose  to  use  wind  magni¬ 
tude  and  direction  as  the  variable  parameters  of  the  problem.  This  approach  is 
followed  in  the  present  study.  The  effects  of  altitude,  Mach  number,  temperature 
distribution,  wind  magnitude  and  direction,  on  sonic  boom  corridor  are  investi¬ 
gated.  #In  the  mathematical  analysis,  a  piecewise  linear  atmospheric  model  is 
assumed  in  order  to  permit  a  closed  form  solution.  The  wind  direction  is  per¬ 
mitted  to  vary  from  altitude  to  altitude.  The  analysis  will  neither  include  the 
effects  of  atmospheric  turbulence  nor  the  effects  of  aircraft  maneuvers.  The 
latter  effects  will  be  taken  up  under  a  separate  study. 

2.0  MATHEMATICAL  THEORY  AMD  ANALYSIS 


The  equations  used  for  calculating  the  propagation  of  sonic  disturbance 
from  a  supersonic  aircraft  are  based  on  the  laws  of  geometrical  acoustics.  These 
laws  state,  in  part,  that  a  wave  front  carrying  a  disturbance  from  a  surface 
(in  this  case,  the  aircraft  body)  of  arbitrary  shape  moves  such  that  its  normal 
velocity  relative  to  the  medium  is  the  undisturbed  speed  of  sound.  The  normals 
are  the  orthogonal  trajectories  of  the  successive  positions  of  the  wave  front, 
and  the  ’rays’  may  be  thought  of  as  carrying  the  disturbance.  In  a  quiescent 
medium,  the  ray  and  the  wave  normal  are  coincident,  but  in  a  moving  medium  (such 
as  atmosphere  with  winds)the  ray  and  the  wave  normal  are  different.  However,  the 
horizontal  component  of  the  wave  normal  velocity  vector  remains  invariant  along 
the  ray. 


In  general,  the  envelope  of  the  acoustic  disturbances  herein  after 
referred  to  as  the  wave  front  (attached  to  the  aircraft  at  the  reference  point 
on  the  aircraft-the  nose)  does  not  have  the  form  of  a  regular  Mach  cone  except 
in  the  special  case  of  straight  flight  at  constant  speed  in  an  atmosphere  of 
uniform  temperature.  The  wave  front  and  che  rays  are  shown  schematically  in 
Fig.  l(a)  for  homogeneous  and  inhomogeneous  atmospheres.  The  calculation  of  the 
geosetry  of  the  rays  as  the  disturbances  propagate  along  them  from  the  aircraft 
to  the  observer  at  the  ground  is  often  called  ’Ray  Tracing’,  and  the  appropriate 
ray  tracing  equations  will  be  presented  in  the  next  section. 

2.1  Acoustic  Ray-Tracing  Equations 

In  order  to  simplify  the  development  of  appropriate  equations  of  ray 
acoustics  (often  called  geometrical  acoustics)  for  our  present  problem,  we  shall 
make  the  following  assumptions: 

(i)  the  speed  of  sound,  c,  and  the  wind  -velocity  vector  W, 
are  independent  of  time  in  the  interval  for  which  ray 
propagation  from  its  origin  to  its  destinatic:.  is  con¬ 
sidered; 

(ii)  wave  front  propagation  occuYs  over  a  sufficiently  small 
region  of  the  earth’s  surface  for  curvature  of  the  earth 
to  be  neglected;  over  this  region  of  space,  the  sound 
speed  and  wind  velocity  vector  are  horizontally  strati¬ 
fied,  i.e.,  no  horizontal  variation  of  sound  speed  and 
wind  velocity  vector  is  permitted. 


The  ray  tracing  equations  for  sound  propagation  from  a  supersonic  air¬ 
craft  sure  conveniently  formulated  in  terms  of  three  coordinate  systems.  The 
coordinate  system  |,i],  z,  positioned  such  that  the  wave  normal  Ires  in  the  l-z 
plane,  is  employed  in  the  analysis  of  ray  geometry;  the  x,  y,  z  coordinate  is 
an  instantaneous  coordinate  for  referencing  the  position  of  the  aircraft  while 
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the  x 


z  is  the  ground-fixed,  observer  coordinate. 


Now  consider  an  aircraft 
z  coordinate  system 


g'  V 

propagating  in  the  negative  direction  of  x-axis  in  the  x,  y, 
and  generating  a  cone  of  disturbance  as  shown  in  Fig.  1(a).  Also  shown  in  Fig. 
i(a)  is  the  ray  cone  whose  tangent  is  inclined  at  the  complement  of  the  Mach 
angle  to  the  aircraft  flight  axis.  Let  us  suppose  that  a  disturbance  was  emitted 
from  the  reference  point  of  the  aircraft  (the  nose)  at  time  =  0,  then  after  a  unit 
tim*,  the  disturbance  would  have  propagated  to  a  terminal  point  along  the  ray  as 
shown  schematically  in  Figs.  l(b),  1(c). 


With  this  understanding,  for  a  particular  ray  under  consideration  the 
propagation  equations  as  derived  in  any  of  Refs.  1  to  k  are  briefly  summarized  as 
follows,  using  the  symbols  in  Figs.  1(a),  1(b),  l(c): 


=  (ic  +  u)  (nc)-1 


g  =  v(nc)'1 


dt 

dz 

ds 

dz 


(nc) 


-1 


1/2 


c/i  +  u  =  ch/jCh  +  u^  =  K,  Snell’s  const. 


(1) 


It  should  be  remarked  that  EG.  (l)  may  be  readily  deduced  from  the 
geometry  of  Figs.  1(b),  1(c),  by  application  of  Snell's  law  for  refraction  of  sound 
at  an  interface  of  discontinuity. 


The  following  brief  discussion  of  Figs.  1(b),  1(c)  is  offered  in  expec¬ 
tation  that  it  may  illuminate  further  understanding  of  the  subsequent  equations. 
In  Fig.  l(b),  the  wind  velocity  vector,  W,  is  shown  in  its  true  magnitude  (in 
accordance  with  the  assumed_horizontally  stratified  medium)  with  its  components 
u,  v.  The  velocity  vector  s^  represents  the  horizontal  projection  of  the  ray 

velocity  vector  s  and  is  displaced  from  the  direction  of  a  unit  normal  vector  n 
by  the  transverse  component  v,  of  wind  velocity.  The  wave  normal  plane  is  in¬ 
clined  to  the  aircraft  flight  axis,  x,  at  an  angle  7  (Fig.  1(b)  )  and  this 
inclination  remains  constant  for  the  ray  under  consideration  throughout  its_  _ 
propagation  time.  In  Fig.  1(c),  Wv,  s^  are  the  projections  of  the  vectors  W,  s  . 
respectively  on  the  %-z  plane  containing  the  wave  normal.  If  the  wind  vector 
is  considered  to  be  directed  along  the  l-axis,  then  W_  =  W.  From  Figs.  l(b), 
l(c)  it  is  readily  seen  that 


I 

v 


s  cos 


c  cos9  +  u 


1(a) 


where 


s  =  cn  +  W 


1(b) 


The  wavefront  is  defined,  according  to  geometric  acoustics,  to  be  propagating  at 
the  undisturbed  speed  of  sound  relative  to  the  medium  through  which  it  passes. 
However,  relative  to  the  ground,  the  velocity  of  the  wave  front  in  its  normal 


\ 


\ 
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direction  is  augmented  by  the  wina^to  yield  (Fig. 1(c)) 


Kc) 


N  =  cn  +  w  .  n 

»  \ 

The  Snell's  constant,  K,  shown  in  Fig.  l(c)  is  readily  seen  to  be  given!  by 

K  =  |n|  sec0 


or 


K  =  (c  +  W  .  n)sec9 


1(d) 


.The  importance  of  Snell’s  constant  in  ray  tracing  calculations  cannot  be  over¬ 


emphasized.  The  ray  of  special  interest  in  the  present  calculations  is  that  ray 
which  reaches  the  ground  at  grazing  incidence  -  the  ray  for  which  the  Snell's 
constant  is  a  maximum  (Sqn.  1(d)),  at  the  ground.  In  the  section  that  follows, 
we  shall  solve  explicitly  the  relevant  ray  acoustic  equations  for  a  chosen  atmos¬ 


pheric  model. 


2.2  Closed  Form  Solution  of  Ray  Acoustic  Equations 


We  observe  from  Fig.  1(b)  that  in  the  presence  of  symmetrical  winds. 


ie.,  head  and  tail  winds  (Wy=  0),  there  is  no  lateral  displacement  of  the  ray  origin 
from  the  wave  normal  plane/  For  oblique  winds  however,  the  ray  is  displaced  from 


, the  wave  normal  plane  in  proportion  to  tne  magnitude  of  the  transverse  component, 
'V,  of  wind  velocity.  For  a  constant  magnitude  oblique  wind  profile  (no  wind 
gradient)  the  effect  of  the  cross  wind  is  merely  the  translation  of  the  wave 
normal  plane  a  distance  v.t  downwind  where  t  is  the  disturbance  propagation  time 
from  flight  altitude  to  the  ground  level.  For  an  oblique  wind  decreasing  to  zero 
velocity  at  ground  level,  the  ray  propagates  down  a  curved  vertical  surface  that 
coincides  with  the  wave  normal  plane  at  the  ground  level  where  v  =  0. 


Fortunately,  the  tangential  component ,  v,  of  wind  velocity,  being  per¬ 
pendicular  to  the  wave  normal  direction  has  no  effect  on  the  application  of 


Snell's  law  to  the  wave  front.  Since  the  Snell's  law  applies  to  the  effective 


wave  normal  velocity  given  by  Ean.  1(c)  and  not  to  the  ray  itself,  ray  tracing 
calculations  are  often  performed  with  the  wave  normal  parameters,  the  ray  para¬ 
meters  being  implicitly  tied  to  it.  Further,  by  suitably  translating  the  wind 
profile  so  that  W  =  0  at  the  ground,  we  may  obtain  ray  and  wave  front  ground  inter¬ 
sections  simply  by  considering  only  the  propagation  in  the  wave  normal  plane  as 
dictated  by  Snell's  law.  For  corridor  calculations,  this  is  exactly  what  is 
done.  Thus  for  the  purpose  of  corridor  calculations  only,  we  make  use  of  the 
first  of  Equations  (l),  the  Snell’s  law,  and  the  direction  cosiue  relation.  , 

Eliminating  the  direction  cosines  £,  a  oetween  the  Snell’s  law  and  the  first 
of  Equations  (l)  gives 


ds 


2.i 


-1% 


(  (k~u)  ~c  )2  (c  +  u(k-u)c  ) 


(2) 


,  The  +  sign  in  Equation  (2)  alxows  for  upward  and  downward  propagating  soniv. 
disturbances.  Equation  (2)  describes  the  space  trajectory  of  a  ronic  disturbance 
propagating  downwards  from  its  source  in  the  ’pave  normal  plane.  The  -  sign 
applies  to  the  downward  propagating  disturbance.  Note  that  the  only  requirement 
for  Equation  (2)  to  correctly  describe  the  ray-ground  intersection  is  that 
W  =  0  at  the  ground  -.  a  situation  which  is  easily  created  by  translation.  v> 
have  used  the  condition  W  =  0  instead  of  v  =  0  (W  does  not  have  to  be  zero  for 
v  to  be  zero)  because  the  whole  analysis  is  carried  out  in  terms  of  wind  magni¬ 
tude  rather  than  its  components.  Thus  the  use  of  W  =  0  rather  than  v  =  0  is 


^  * 


for  symbolic  convenience  rather  than  a  limitation.  It  must  be  pointed  out  that 


for  consideration  of  the  geometry  of  the  ray,  through  the  atmosphere ,  ray  propa¬ 
gation  in  ^-direction  must  be  taken  into  account.  For  conpleteness  of  this 
study,  the  expression  for  ray  propagation  in  ^-direction  is  derived  herein  at 
the  appropriate  point. 

1  If  the  .speed  of  sound  c  and  wind  velocity  component  u  along  the  wave 
normal  plane  are  known  as  functions  of  z,  then  Equation  (2)  may  be  evaluated 
directly  by  quadrature.  However,  it  is  instructive  to  assume  a  piecevise  linear 
atmosphere  and  seek  a  closed  form  solution  of  Equation  (2).  Accordingly,  define 
the  model  wind  profile  ■  * 

u'  =  ug  +  az 

and  sound  speed  profile 

c  =  eg  -  pz 


so  that  relative  to  the  wind  speed  at  the  ground  ug,  we  have 

i  u  =  az 

c  =  eg  -  pz 


(3) 


where  a,  p  are  positive  quasi-constants  (piecewise  constant  wind  and  sound  speed 
gradients  respectively).  Upon  substituting  Equation  (3)  into  Equation  (2)  we 


obtain  the  ray  slope  equation  as 


(«0 


Let 


Then 


where 


z^  =  k  -  az 

z2  =  Ypz 

zr  ’  Vz2 

dz 
d£ 


ac  z  -  akz  +  s 
g  r  r  o 


(5) 


(6) 


s  =  pk  -  ac 
°  g 

Combining  the  paramet'_rs  in  Equation  (5)  and  differentiating,  gives: 


dz  = 


s  dz  ‘ 
o  r 


(pzr-  a)2 


and  using  this  in  Equation  (6)  gives: 


5 


(7) 


dz  [z  2-l]2(pz  -a)2 

_£  =  _  _£ _ _  r  7  ,  (7) 

d|  2  ,  'l} 

ccc  z  -akz  4  s 
g  r  r  o 

The  steps  leading  to  Equations  (5)  -  (7)  are  necessary  to  bring  Equation  (^4)  to  a 
form  permitting  the  use  of  standard  integrals.  Intergrating  Equation  (7)  we  have 


6“  '0  = 


ac  z  -  akz  +  s 
g  r  r  o 

Ozr-a)2(zr2-i)1/2 


which  upon  employing  partial  fractions  integrates  to 

6  -  60  =  -  {(Pk-ocg)l4+  a/P[ai4+  -  k)  } 


where 


1  1  r  (P2-a2)-a(Pz^-a) 

I-  =  — ~~ p-r  tan  x - - - 

3  (aW  !  P'j(a2-p2)(z  2.1}] 


1  r  e(*r-1)2~ 

4  =  -  p--p  ai  -  - £ - 

(P  -a  )  3  (Pzr-a) 


k  -  az 

2r  c  -Pz 
g 


After  considerable  algebraic  manipulations,  Equation  (9)  reduces  to 

5i*i  '  5i  +  {  f(k-Vi)2-c?n] 3  -  [(t-ui>2-  4  j  3  } 


x  (z.-z.  .)  (c  -c.) 
'  1  1+1  1+1  1 


Strictly  speaking.  Equation  (13)  represents  the  true  ray  trajectory  in  the  wave- 
normal  plane  only  for  symmetrical  winds.  For  oblique  winds  decreasing  to  zero 
velocity  at  the  ground,  Equation  (13)  represents  the  projection  of  ray  trajec¬ 
tory  on  the  wave  normal  plane  until  the  ground  is  reached  at  which  point  it  repre¬ 
sents  the  true  ray  trajectory. 

Noting  that  Equation  (13)  becomes  singular  when  c^+^  =  c^,  we  seek  a 

solution  valid  under  this  condition.  Physically,  this  condition  inrolies  that 
sonic  boom  can  spread  to  infinite  distances  in  an  isothermal  atmosphere.  The 
solution  valid  under  isothermal  conditions  in  the  nresence  of  winds  is 
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I  = 


i  I  {  44'Vi1  +  f  (k-ui+1)2-4i|2  |  ♦  k 

i=l 

(k-ut)  +  [(k-u.)2  -  c2nj  2 


%2  2  * 
(k-Vl)  -  oi+1^ 


-  ci+l  in 


(k-Ui)2-  c2+1j  2  } 


X  {z.-  Z.)  c7^,  (u.,  -  u.) 
'  1+1  l'  1+1  1+1  l' 


-1 


m 


N  =  number  of  atmospheric  layers  for  which  data  is  given. 


Under  isothermal  conditions  and  no  winds,  the  valid  solution  derived  from  Equation 
(13)  is: 


|  =  c  (z.-z  }  [k2-  c2  f 1/2 

*  l+l  v  l  l+l*  1  i+lJ 


(15) 


As  discussed  in  Section  2.1,  for  symmetrical  winds  the  ray  is  contained  in  the 
wave  normal  plane,  hence  the  equations  developed  up  to  this  point  are  sufficient 
to  carry  out  ray  tracing  and  corridor  width  calculations.  However,  for  complete¬ 
ness  of  this  study,  and  to  enable  the  reader  perform  ray  tracing  in  the  presence 
of  cross-winds,  the  ray  propagation  in  ^-direction  must  be  included.  To  do  this, 
we  eliminate  the  direction  cosines  £,  n  between  the  second  of  Equations  (l)  and 
the  Snell's  law,  to  obtain 

=  _  vc-1  (k-u)  [(k-u)2  -  c2]  ~1,/2  (16) 


The  velocity  components  u,  v  are  related  (see  Fig.  1(b))  through  the  expressions 
u  =  -  Wcos(H-D-?) 

v  =  Wsin(H-D-?)  (17) 

so  that 

v  =  -u  tan(H-D-7)  (18) 

Substituting  Equation  (l8)  in  (16)  and  going  through  the  same  procedures  as  before, 
one  obtains 


= 

dz 


Ct  z  tan(H-D-7) 
(cg-Pz) 


(k-az) _ 

[(k-u)2-(cg-pz)cj  1/f2 


which  after  substitutions  and  use  of  partial  fractions  integrates  to 


(19) 
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v1”1  rr  -i  rci+r^^K'ui+i^ 

*  I  “iKV*)  {  [tni+l-T“  f[(K.u  )2.02  ]l/2 

i=l  LV  ui+l;  i+lJ 


-  ^ 


r/ir  \2  2  ,1/2 
(K"Ui+l“4iCi+l) 


r  -i  rci"ui(K'ui)  i  L  lK~  v  "ci J '  ii 

- 1 *v™  777  ,2  W ■■  -  —  } 

l[(H)  -Cj]  '  J  (K-Ui-Jiici) 


r/l.  \2  2  ,1/2 
[(K-^)  -ci  ]  • 


x  (Z.,,-  Z.).(c  -  c.) 

v  l+l  l  l+l  i7 


-1 


(20) 


where 


"i  *  (Vl  -  V  <ci+l  -  ci> 


-1 


tn.+1  =  tan  (H-D.+1  -  y) 


(21) 


Under  isothermal  conditions,  p.  aqd  Equation  (20),  like  Equation  (13) >  is 

invalid.  The  valid  solution  is  then 

N-l 

’l  =  "  -1/2  tni+l  ^k+  ui+i)[(k~  ui+i^  "  ci+i^  ^  "  ci+l^n  ^k~  Ui+1^ 

i=l 

+  (k-  u1+1)2-  c2+1  ]1/2|  J  '  tDi  {<>'-"!>[  (k-“i>2  -  4  11/2 
-  4n  ini<k-ui) +  [<k-v2  -  4+i]1/2  ] } 


*  <zi-l  -  <ui+l  -  ui> 


-1 


(22) 


The  lateral  spread  of  sonic  boom  is  obtained  by  coordinate  transformation  - 


x  =  |cos7  -  qsin7 

y  =  |sin7  +  T)cos7  (23) 

The  maximum  lateral  spread  of  the  sonic  boom  is  given  by  the  so  called  ’cut-off’ 
point.  This  point  is  determined  by  the  last  ray  to  strike  the  ground  at  grazing 
incidence.  The  necessary  condition  and  equation  for  fulfilment  of  grazing  inci¬ 
dence  will  now  be  considered. 
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2.3  Geometrical  Considerations 


The  purpose  of  this  section  is  to  express  the  wind  components  along  and 
perpendicular  to  the  wave  normal  plane  in  terms  of  the  wind  magnitude  and  direc¬ 
tion,  and  to  define  the  wave  normal  plane  for  the  condition  of  grazing  incidence 
for  the  ray  under  consideration.  The  angle  defining  the  wave  normal  plane  is 
needed  to  transform  the  wave  normal  coordinate  system  to  the  aircraft  coordinate 
system. 

Consider  the  aircraft  to  be  moving  in  the  negative  direction  of  x-axis 
as  shown  in  Pig.  l(a).  Then,  relative  to  the  coordinate  system  moving  with  the 
wind  at  the  aircraft  altitude,  the  following  geometrical  relations  are  derived 
from  Fig.  1(a): 

sin7  =  Asin0  [1  +  A2sin2  0]  ^2 


.cose  =  -M"1  [1  +  ?v2sin20  ]-1/2 


where 


cos0  =  -(Mcos7) 
X-  (M2-!)'1/2 


It  is  convenient  to  perform  the  calculations  in  a  coordinate  system  moving  with  the 
wind  because  it  lends  itself  to  the  use  of  simple  geometrical  relations  such  as 
given  in  Equations  (24).  If  calculations  are  performed  in  the  coordinate  system 
fixed  to  the  earth,  then  the  deformation  of  the  initial  wave  system  due  to  the  wind 
at  the  flight  altitude  will  have  to  be  considered.  In  the  coordinate  systan 
moving  with  the  wind  at  the  flight  altitude,  the  wind  velocity  is  zero  at  the 
flight  altitude.  The  final  result  required  in  a  ground  fixed  coordinate  system 
is  obtained  through  a  galilean  transformation. 

We  shall  specify  the  wind  direction  such  that  the  heading  angle,  H,  is 
measured  clockwise  positive  from  the  true  Worth  and  winds  reported  as  blowing  from 
the  direction  D,  also  measured  positive  clockwise  from  the  Worth  -  Fig.  1(b). 


path  are: 


The  components  of  wind  velocity  along  and  perpendicular  to  the  flight 


Wx  =  -Wcos(H-D) 


Wy  =  Wsin(H-D) 

hence  the  component  of  wind  velocity  parallel  to  wave  normal  plane  is 

u(z,7)  =  Wx  cos7  +  Wysin7 
=  -Wcos(H-D-7) 
and  perpendicular  to  wave  normal  plane  is 

v(z ,7)  =  -Wx  sin7  +  Wy  cos7 
=  Wsin(H-D-7) 
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The  angle,  7,  defining  the  particular  wave  normal  plane  for  grazing  incidence 
at  the  ground  may  now  be  determined.  Rewrite  Snell’s  law  given  in  Equation  (l) 
as: 


K'  =  c  /cos0  +  u,  =  cg/cose  +  ug 
B  B  B  S 


(28) 


The  condition  for  grazing  incidence  at  the  ground  requires  that  cos0  =  -1, 
so  that  from  Equation  (28),  we  have  ® 

K'  =  -eg  +  ug 

Then  relative  to  the  coordinate  system  moving  with  the  wind  at  the  aircraft 
altitude,  we  have  -  C, 


and 


k  =  --eg  -  u^  +  ug 
C°s0h  =  Cjj/k 


(29) 

(30) 

(3D 


Note  that  by  virtue  of  the  use  of  a  moving  coordinate  system,  the  ray  and  the 
wave  normal  are  coincident  at  the  flight  level  where  n  =  0,  hence  the  use  of 
initial  wave  normal  slope  in  this  calculation.  From  geometrical  considerations 
(Fig.  1(b)) 

u(z,  7)  =  -Vicos (H-D-7) 


hence 


cos9h  =  cfa[-cg  +  (W  -Wh) cos (H-D-7)] 


-1 


(32) 


But  from  Equation  (24) ,  cosB^  =  -(M  COS7)'1  and  using  this  in  Equation  (32)  gives 


2,  t-2 


2  4 

ip 


c  (c  M-Wcosw)  +[Wsinu  (c  M)  +  w  -  2c,  MWcosu  -  c  ]2 


cos7  = 


g  n 


(c,M)2  -  2c,  Mfcosw  +  W2 
v  h  '  n 


fL 


(33) 


where  u  =  H-D.  The  -  sign  before  the  radical  applies  to  the  leeward  side  and 
the  +  sign  to  the  windward  side  of  the  flight  track. 


The  angle,  7,  defining  the  wave  normal  plane  is  given  by  Equation 
(33)  while  the  initial  wave  normal  slope  at  emission  altitude  is  given  by 
Equation  (32).  Equation  (33)  is  important  and  leads  to  some  interesting  results. 
It  is  possible  to  determine  the  magnitude  of  wind  velocity  required  to  confine 
sonic  boom  to  the  flight  path  alone, 

2.4  Winds  Required  for  ’Cut-Off1  of  Sonic  Boom 

Under  the  flight  path,  7=0,  0  =  n’,  the  conrolement  of  the  Mach  angle 
relative  to  wind.  If  the  heading  angle  is  taken  arbitrarily  as  l80°  (Fig.  1(b)) 
then  for  head  wind,  D  =  l80°  so  that  H-D  =  0.  Substituting  these  parameters  in 
Equation  (33)  gives  the  limiting  head  wind  as: 

W  =  chM  -  eg  (34) 
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Similarly  for  tail  wind. 


2.5  Ray  Propagation  Time  Histor 


The  ray  propagation  time  is  given  in  Equation  (l)  as: 
dt  ,  ,-l 


where 


2  2-,  k  ,-l 


nc  =  c[(k-u)  -c  ]  2  (k-u) 


Substituting  Equation  (37)  in  (36)  and  integrating  between  any  two  points,  gives 
the  ray  propagation  time  between  these  two  points  as: 


t  =  f  (k-u)  [(k-u)2-c2]  _a  c-1 


Making  use  of  Equation  (3)  in  (38)  gives: 


t  =  f  (k~az)(c  -3z)_1[(k-OJz)--(cg-Bz)2  ]~2  dz 


If  we  go  through  the  same  arguments  and  substitutions  as  in  Section  2.1,  we  obtain 
from  Equation  (39)  >  the  following  result: 


T  r  _i  p  I 

th-t  =  I  Zr(Szr-a)  l[zrd-l]  -2 


where  is  the  time  of  emission  of  the  ray  at  the  flight  altitude,  and  z  is 
defined  in  Equation  (5). 


By  writing 


**  B  B'1  [*r2-lP  {l+  a/p^-a/e)'1} 


and  substituting  in  Equation  (^0),  one  obtains 


-l  -1  -?r  ? 

t^-t  =  B  cosh  zp+  ap  I  [l-(a/B)  ]  2TAN  j* - 


-(a/3)  )-a/P(z _-a/P) 


tl-(a/e)2)(2r2-l)l 


m  1  -1 1  r 
=*">]]  !H) 


Let  q  =  a/B  as  defined  in  Equation  (21),  then  Equation  (itl)  becomes 


V‘  -  ^  f-*'1  (¥  )  *  -  ™  1  [fg^TTB  ]  } 


where  it  is  understood  that  in  Equation  (*42)  u  =  u(z),  c  =  c(z).  Recasting  Equation 
(42)  into  a  form  more  suited  to  digital  operations,  one  obtains: 

t .  y  {cosh-1  ( i2sa)  ♦  | 

tl  1  V  i^  '  1  l[(K-u.+1)2-  o\^/2  1 

_ “i  (  K-Ui\  ..  „.„-i  rv'i^-V  1 1 


v  \j  1  l[(K-u/-  c2  ]V2  IS 


Note  that  Equation  (43)  represents  a  gener  :1  equation  for  calculation  of  propa¬ 
gation  time  for  any  ray  which  may  or  may  nt  reach  the  ground.  For  a  quiescent 
state  between  layers.  Equation  (43)  reduces  to 

N-l 

t  =  £  {cosh*1  -S-  -  cosh'1  m 

i=l  1 

For  a  special  case  in  which  a  ray  is  required  to  reach  the  ground  at  grazing  inci¬ 
dence,  Equation  (44)  takes  the  form 

t  =  P  1  cosh-1  (cg/ch)  (4f/) 

where  P  is  the  sound  speed  gradient. 

This  shows  that  for  this  particular  ray,  the  propagation  time  from  the 
flight  altitude  to  the  ground  depends  only  on  the  altitude,  for  a  given  ground 
condition. 

The  general  solution  applicable  to  isothermal  state  between  layers  with 
winds  is  given  by 

1  (ci«“irl{  [k-\)2  -  4-  4i]1/2-i<k-ui+/-  4ii1/2} 


For  an  isothermal,  quiescent  state  between  layers,  the  ray  propagation  time  is 
simply 

t  =  k  c-y  (h2  -  c2+1)-V2  (  Vl-  z.)  (XT) 

The  transformations  between  a  ground  fixed  coordinate  system  and  that  defined  by 
Equation  (23)  are 
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x  =  X  -  W.cos(H-D-.)  t 

8  h  g  (1,8) 

yg  =  ¥  +  Whsin  (H-D-  -)t 

2.6  Shock-Ground  Intersection 

Equation  (^8)  gives  the  raywground  intersection  when  z  =  0.  We  must 
clearly  distinguish  the  difference  between  the  ray  and  the  shock-ground  inter¬ 
sections.  By  ray-ground  intersection,  we  mean  the  ground  locus  of  acoustic  dis¬ 
turbances  emitted  from  the  aircraft  simultaneously.  These  disturbances,  how¬ 
ever,  arrive  at  the  ground  at  different  times,  t  ,  depending  on  their  initial 
direction  of  propagation  at  emission, time,  tu.  ®  The  shock-ground  intersection 
on  the  other  hand  is  the  locus  of  sonic  disturbances  reaching  the  ground  simul¬ 
taneously  -  as  sonic  boom.  These  disturbances  may  have  been  emittbd  at  different 
times. 

For  steady,  level  flight  being  treated  herein,  the  coordinates  of  shock 
-ground  intersection  are 

X  =  x  —  W  cos(H-D-  )t  — c  M(t  -t  )  (h9) 

s  hv  'mh'gnr  K 

y  =  y  +  W,  sin(H-D  -  ;)t 
*  h  v  ' '  m 

where  t  =  time  of  arrival  of  disturbance  at  the  ground  along  the  ray 

o 

tffl  =  time  of  arrival  at  the  ground,  of  the  disturbance  propagating 
directly  under  the  flight  path  (minimum  ray  propagation  time). 

We  shall  in  the  next  section  discuss  the  results  obtained  from  computer  solu- 
t-'ons  of  the  equations  developed  in  this  section. 

3.0  DISCUSSION  OF  RESULTS 

3.1  Standard  Atmosphere 


In  the  last  section  we  discussed  the  significance  of  the  analytical 
expressions  on  a  qualitative  basis.  In  this  section,  we  shall  discuss  the 
computational  results  on  a  qualitative  basis.  It  is  reasonable  to  base  our 
discussions  on  flight  altitudes  within  the  stratosphere,  since  the  proposed 
flight  altitudes  of  the  SET’s  fall  within  this  domain.  Now,  for  a  given 
atmospheric  model,  the  main  parameters  among  others,  that  influence  the  sonic 
boom  corridor  are  the  flight  altitude  and  Mach  number,  wind  and  temperature 
gradients  and  wind  direction.  It  is  therefore  natural  to  determine  the  varia¬ 
bility  of  sonic  boom  corridor  in  terms  of  these  physical,  measurable  parameters. 

Equations  (3*0,  (35),  displayed  graphically  in  Fig.  2  show  that  for 
flights  above  the  tropopause  in  a  standard,  quiescent  atmosphere,  sonic  booms 
do  not  reach  the  ground  for  M  =  C  /C,  <  1.15.  This  Mach  number  apparently 
corresponds  to  the  so-called  ’cut-ofr'  Mach  n’unber  in  a  quiescent  atmosphere. 
The  term  ’so-called'  is  used  here  because  booms  are  still  heard  on  the  ground 
until  well  past  the  ’cut-off’  condition.  At  this  condition  for  which  the  wave 
cusp  just  reaches  the  ground  one  has  a  superboom  (focussed  boom)  on  the  ground. 
As  we  move  further  into  the  'cut-off  region,  the  cusp  moves  higher  above  the 
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ground  and  the  boom  degenerates  into  a  mere  rumble.  Beyond  the  ’cut-off’  point, 
the  effect  of  the  acoustic  disturbances  is  negligible. 

A  further  interesting  characteristic  of  the  tail  wind  is  clearly 
evident  in  Fig.  2.  Observe  that  whereas  ’cut-off’  Mach  number  for  flights 
above  the  tropopause  in  a  quiescent  atmosphere  is  1.15,  tail  winds  ranging  from 
zero  to  87.6  knots  at  the  flight  level  can  reduce  the  ’cut-off'  Mach  number  from 
1.15  to  1.0.  F\irt. her  remarks  regarding  the  interpretation  of  Fig.  2  will  be 
made  when  discussing  the  effects  of  altitudes  on  sonic  boom  corridor. 

The  effects  of  flight  Mach  number  and  wind  speed  on  sonic  boom  corri¬ 
dors  are  shown  in  Fig.  3.  In  this  figure,  five  groups  of  corridors  (length  of 
each  strip  represents  corridor  widths;  width  of  each  strip  is  non-dimensional) 
labelled.  A,  B,  C,  D,  E,  are  shown.  Within  each  group,  starting  from  the  left 
and  proceeding  to  the  right,  the  strips  represent  no  wind,  tail  wind,  side  wind 
(blowing  from  right  of  flight  path  to  the  left) ,  and  head  wind  respectively. 

The  coordinate  position  of  each  group  of  corridors  in  relation  to  the  zero  corri¬ 
dor  lines  -  PQ,  PR,  is  similarly  labelled  in  Fig.  2,  in  a  wind-Mach  number  cor 
ordinate  system.  Since  line  PR,  Fig.  2,  represents  wind  and  Mach  number  com¬ 
binations  for  which  the  corridor  is  reduced  to  a  line  along  the  flight  track 
for  head  winds,  it  is  to  be  expected  that  as  we  move  away  from  IR  along  a  line 
of  constant  wind  (as  points  A,  B,  C)  the  corridor  width  should  increase,  as 
indeed  it  does  in  Fig. 3-  An  alternative  description  of  the  line  PR  in  Fig.  2 
is  that  it  represents  wind  speed  and  aircraft  Mach  number  combination  for  which 
only  the  ray  directly  under  the  flight  path  just  reaches  the  ground  tangentially. 
Thus  for  an  aircraft  in  a  straight,  level  flight  at  constant  speed,  the  rays 
emitted  at  successive  times  will  form  an  envelope.  PR  and  IQ  are  thus  superboom 
paths  for  head  and  tail  winds  respectively  under  the  aircraft. 

The  author  shares  the  opinions  expressed  in  Refs.  3>  5  that  the  ratio 
’corridor  with  winds/corridor  without  winds’  remains  essentially  constant  for 
aircraft  speeds  larger  than  M  =  1.3.  However,  the  author  wishes  to  stress  that 
the  opinions  expressed  in  ‘-.hese  references  -  that  the  ratio  ’corridor  with 
winds/corridor  without  winds*  is  close  to  unity  for  M  >  1.3  is  valid  only  for 
small  and  medium  wind  profiles  as  used  in  those  references  and  shown  in  Fig.  10. 

It  happens  that  the  magnitude  of  the  wind  speed  at  the  flight  altitude 
has  a  predominating  influence  on  the  corridor  width,  over  the  effects  of  actual 
wind  distribution,  just  as  the  effect  ground  temperature  predominates  that  of 
temperature  distribution.  Now,  for  the  'mean  zonal'  wind  profile  used  in  Ref.  3> 
at  50,000  ft.,  the  wind  speed  is  about  31.4  knots  and  at  the  ground  is  5.9  knots. 
Thus,  relative  to  wind  at  the  ground  fixed  coordinate,  the  wind  at  flight  alti¬ 
tude  is  25.5  knots,  which  is  small  compared  with  the  speed  of  sound  at  the  flight 
altitude .  In  the  presence  of  strong  winds,  such  as  200  knots,  the  ratio  'corridor 
with  winds/corridor  without  winds'  is  neither  constant  nor  close  to  unity  as 
is  shown  by  corridor  groups  D  and  E  in  Fig.  3» 

Further  effects  of  wind  magnitude  and  direction  are  shown  in  Fig.  4 
for  constant  altitude  -  40,000  ft.  and  M  =  1.6  for  linear  wind  profiles  decreasing 
to  zero  at  the  ground.  Observe  that  for  a  20-knot  wind,  the  increase  in  corridor 
width  by  a  tail,  wind  (increase  above  no  wind  corridor)  is  about  the  same  as  the 
decrease  by  a  head  wind  -  the  type  of  effect  expected  from  a  linear  theory  solu¬ 
tion.  However,  as  the  wind  speed  increases,  the  departure  from  the  above  trend 
becomes  more  pronounced.  The  effect  of  side  (cross)  winds  (D  =  90°)  is  clearly 
evident  in  Fig.  4.  The  general  effect  of  the  side  winds  is  the  displacement  of 
the  boom  corridor  to  the  lee  of  the  flight  track.  The  particular  effect  is  a 


slight  reduction  in  corridor  width  for  strong  winds  (W  >  60  knots,  say).  Observe 
in  Fig.  4  that  for  a  side  wind  of  20  knots,  the  corridor  width  is  unaltered  from 
the  no  wind  case.  For  W  =  64.5  knots,  the  corridor  width  is  reduced  by  2  miles  ■ 

from  the  no  wind  case,  and  for  W  =  200  knots,  by  5  miles  or  13$.  Thus  for  small  j 

and  moderate  side  wind  profiles,  the  boom  corridor  width  is,  for  all  practical  j 

purposes,  the  same  as  for  no  wind.  \ 

The  effect  of  flight  altitudes  on  sonic  boom  corridor  is  calculated  and  j 

shown  in  Fig.  5  for  M  =  1.6,  W  =64.5  knots  at  the  flight  altitude  and  decreasing  j 

linearly  to  zero  at  the  ground.  This  figure,  which  is  self-explanatory,  in  effect  \ 

illustrates  the  effect  of  wind  gradients.  High  wind  gradient  reduces  boom  corri-  \ 

dor  widths,  just  as  high  temperature  gradient  does.  We  can  infer  from  Fig.  5 
that  if  a  64.5-knot  wind  at  40,000  ft.  creates  a  narrower  corridor  on  the  ground  ■ 

than  the  same  wind  at  50,000  ft.,  then  a  wind  that  causes  a  ray  emanating  from  j 

50,000  ft.  to  reach  the  ground  only  under  the  flight  track  at  grazing  incidence  j 

(a  point  on  FR,  Fig.  2)  will  cause  complete  cutoff  for  the  ray  emanating  from  i 

4o,OOQ  ft.  With  this  explanation  in  mind,  let  us  now  go  back  to  Fig.  2.  The  ; 

inference  made  from  Fig.  5  implies  that  although  Fig.  2  is  valid  for  flights  j 

between  36,080  ft.  and  65,000  ft.,  the  actual  altitude  of  origin  of  the  distur-  I 

bance  must  be  known  before  a  valid  interpretation  of  Fig.  2  can  be  made  for  any  \ 

case  of  interest.  This  clarification  has  not  been  made  in  any  literature  where 
Fig.  2  has  been  presented,  and  the  author  feels  that  the  absence  of  this  clarifi¬ 
cation  could  lead  the  reader  to  wrong  conclusions. 

The  present  calculation  is  compared  with  the  results  of  Ref.  3  for  the 
mean  zonal  wind  profile  shown  in  Fig.  10  in  a  standard  atmosphere,  at  50,000  ft., 

M  =  1.2.  Agreement  of  both  results  is  excellent,  despite  the  fact  that  Ref.  3 
did  not  assume  linear  wind  and  temperature  profiles. 

3.2  Non-Standard  Atmospheres 

In  general,  ray  propagation  characteristics  depend  on  the  local  varia¬ 
tions  of  atmospheric  properties;  however,  whether  a  ray  gets  to  the  ground  or  not, 
depends  mainly  on  the  distributions  of  wind  and  temperature  at  the  flight 
altitude  and  the  ground.  It  is  inferred  from  Equation  2(d)  that  the  wave  normal 
which  determines  the  maximum  lateral  spread  of  sonic  boom  is  the  one  for  which 
the  Snell's  constant  is  a  maximum.  Relative  to  the  coordinate  system  moving 
■with  the  wind  at  the  flight  altitude.  Equation  (32)  indicates  that  the  initial 
inclination  to  the  horizontal  of  the  wave  normal  that  is  to  reach  the  ground  at 
grazing  incidence,  is  dictated  solely  by  the  wind  and  the  sound  speeds  at  the 
flight  level  and  the  ground.  Thu*'  for  large  temperature  difference,  AT,  between 
the  ground  and  the  flight  altitude,  the  ray  leaves  its  source  at  a  steep  inclina¬ 
tion  to  the  horizontal  and  hence  can  withstand  temperature  inversions  typical  of 
non-standard  atiaospheres .  On  the  other  hand,  for  sms, 11  AT,  the  ray  leaves  its 
source  at  a  shallow  inclination  to  the  horizontal,  thus  making  it  very  sensitive 
to  temperature  inversions.  Under  this  condition,  a  slight  abrupt  change  in 
atmospheric  properties,  such  as  temperature  inversions,  results  in  very  large 
values  for  sonic  boom  corridors. 

For  the  purpose  of  comparison,  sonic  boom  corridors  have  been  com¬ 
puted  for  a  typical  Canadian  town,  Maniwaki,  Quebec,  for  the  month  of  January, 
using  the  mean  meteorological  data  provided  in  Refs.  9  and  10  and  plotted  in 
Figs.  8  and  9«  Comparison  is  made  in  Fig.  7,  with  the  corridors  labelled  group 
C  in  Fig.  3.  Observe  that  with  AT  (winter)  just  about  60$  of  AT  (standard  atm.), 
winter  corridor  is  more  than  triple  the  standard  atmosphere  corridor  for  tail 
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wind,  other  meteorological  parameters  held  constant.  It  is  interesting  to  note 
that  the  temperature  at  the  flight  altitude  -  50,000  ft,  remains  about  the  same 
for  winter  and  standard  atmospheres  (see  Fig.  8),  so  that  the  ground  temperature 
is  mainly  responsible  for  the  difference  in  corridors.  This  leads  us  to  conclude 
that  the  ground  surface  temperature  must  be  the  most  important  meteorological 
factor  affecting  the  sonic  boom  corridor. 

4.0  CONCLUSIONS 

Calculation  of  sonic  boom  corridors  based  on  closed  form  solutions  of 
the  ray  acoustic  equations  using  piecewise  linear  atmospheric  models  of  winds 
and  temperatures  has  been  accomplished.  The  computer  program  for  these  calcula¬ 
tions  processes  5  different  cases,  each  of  8  atmospheric  layers  in  3  minutes 
using  the  comparatively  slow  computer  -  IBM  1130. 

The  results  show  that  whereas  a  complete  cut-off  of  sonic  boom  for 
flights  above  the  tropopause  in  a  quiescent  standard  atmosphere  occurs  for 
M  <  1.15,  the  cut-off  Mach  number  is  reduced  by  the  tail  wind.  Specifically, 
a  tail  wind  of  about  87  knots  will  reduce  the  cut-off  Mach  number  to  M  =  1.0 
at  the  appropriate  altitude  within  the  tropopause. 

Theoretically,  sufficiently  high  head  winds  will  confine  the  sonic 
to  the  flight  track  only,  with  attendant  focussing  effects;  higher  winds  will 
cause  boom  cut-off.  However,  the  required  head  winds  increase  with  Mach 
number  and  exceed  200  knots  above  M  =  1.5.  The  main  effect  of  side  winds 
is  to  shift  the  corridor  laterally  leeward  with  respect  to  the  flight  track, 
the  shift  being  in  preportion  to  the  wind  strength.  Side  winds  less  than  20 
knots  at  40,000  ft.  do  not  alter  the  width  of  the  corridor  from  the  no  wind 
case,  but  higher  winds  cause  a  slight  reduction.  In  particular,  a  side  wind 
of  200  knots  at  40,000  ft.  causes  a  13$  reduction  in  (corridor  width  from  no 
wind  case. 

It  is  found  that  for  small  to  moderate  wind  profiles  and  M  >  1.5, 
the  increase  in  corridor  width  (above  no  wind  case)  due  to  tail  wind  is 
approximately  equal  to  the  decrease  due  to  head  wind.  For  stronger  winds 
at  the  same  Mach  number,  head  winds  produce  progressively  higher  decrease  in 
corridor  width  than  the  increase  due  to  tail  winds.  The  largest  variations 
in  corridor  widths  due  to  tail  winds  occur  for  M  <  1.3. 

The  effect  of  winds  on  sonic  boom  corridor  is  more  pronounced  for 
flights  above  the  tropopause  wnere  isothermal  conditions  prevail,  but  is  less 
significant  for  flights  below  the  tropopause  where  temperature  effects  are 
dominant.  Based  on  the  results  for  non-standard  atmospheres,  the  ground 
temperature  is  the  greatest  single  meteorological  parameter  affecting  the 
sonic  boom  corridor;  the  influence  of  ground  temper at ure  is  such  that  higher 
than  standard  temperature  constricts  it  whilst  lower  than  standard  temperature 
expands  it. 


16 


REFERENCES 


1.  Keller,  J.  B. 
2*  Groves,  G.  V. 


3.  Hayes,  W.  D. 
Haefeli,  R.  C. 
Kulsrud,  H.E. 

4.  Friedman,  M.  P. 


5.  Kane,  E-.  J. 
Palmer,  T.  Y. 

6.  Randall,  D.  G. 


7.  Dressier,  R. 
Fredholm,  N. 

8.  Reed,  J.  W. 
Adams,  K,  G. 


9.  Titus,  R.  L. 


Geometrical  Acoustics  I.  The  Theory  of  Weak  Shock 
Waves.  Jour.  Appl.  Phys.,  Vol.25,  No. 3,  1954. 

Geometrical  Theory  of  Sound  Propagation  in  the 
Atmosphere.  Jour,  of  Atmospheric  and  Terrestrial 
Fhys.  Vol.7,  1955. 

Sonic  Boom  Propagation  in  e.  Stratified  Atmosphere, 
With  Computer  Program.  NASA  CR-1299,  April  1969. 


A  Description  of  a  Computer  Program  for  the  Study 
of  Atmospheric  Effects  on  Sonic  Boom.  NASA  CR-157 
1965/ 

Meteorological  Aspects  of  the  Sonic  Boom.  The 
Boeing  Company,  Renton,  Washington,  September  1964 

Methods  for  Estimating  Distributions  and  Inten¬ 
sities  of  Sonic  Bangs.  R.A.E.  R  &  M  No.  3U3, 
August  1957. 

Statistical  Magnifications  cf  Sonic  Booms  by 
the  Atmosphere.  FFA  Report  No.  104,  June  1966. 

Sonic  Boom  Waves  -  Calculation  of  Atmospheric 
Refraction,  Aerospace  Engrg.  Aerodynamics  and 
Fluid  Mechanics.  March  1962. 

Upper  Air  Climate  of  Canada  -  Average,  Extreme 
and  Standard  Deviation  Values.  1951  -  i960. 

Dept,  of  Transport,  Meteorological  Branch, 

1965. 


10.  Henry,  T.  J.  G 


Maps  of  Upper  Winds  Over  Canada.  Dept,  of 
Transport,  Meteorological  Branch,  1957. 


APPENDIX 


Snell's  Law  and  Ray  Traci] 


The  purpose  of  this  appendix  is  to  define  Snell's  law  in  the  context 
used  for  geometric  acoustic  propagation  in  a  stratified  atmosphere,  and  to 
eiq>lain  its  role  in  ray  tracing.  ! 

In  geometric  acoustics,  disturbances  are  propagated  on  the  wave  front 
along  the  rays.  The  wave  front  propagates  such  that  its  normal  velocity  rela¬ 
tive  to  the  medium  is  the  undisturbed  speed  of  sound.  I 

\  \ 

i  In  geometric  optics,  Snell's  law  of  refraction  defines  the  refractive 
index  of  a. medium  in  terms  of  the  speed  of  propagation  of  light  wave  through 
it,  i.e.. 


Refractive  index  = 


speed  of  light  in  vacuum 
speed  of  light  in  medium 


For  the  propagation  of  a  sound  wave  in  a  horizontally  stratified, 
quiescent  atmosphere,  an  analogous  refraction  law  states  that ithe  trace  velocity 
of  the  wave  front  along  the  interface  separating  two  layers  is  conserved.  Even 
when! there  is  relative  motion  between  layers  (see  sketch  below)  the  conservation 
of  the  trace  velocity  of  the  wave  front  still  holds;  thus  the  equivalent  of 
Snell's  law  for  a  stratfied  moving  medium  may  be  states  mathematically  as 


+  U1  = 


+  u2  =  constant,  k 


A1 


Incident 


not  the'  actual  ray  path  -  a  kinematic  entity,  which  is  used  for  ray 


tracing,  hut  Its  geometric  equivalent  -  its  projection  on  the  wave- 
normal  plane.  Hence  to  calculate  the  ray,  the  information  provided 
hy  Snell’s  law  about  the  wave  normal  slope  at  each  interface  is  used 
to  evaluate  the  slope ^ of  the  geometric  ray  at  the  interface. 
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FIG.  1(a).  SKETCH  SHOWING  COORDINATE  ORIENTATION  AND 
INITIAL  WAVE  NORMAL  DIRECTION  RELATIVE  TO 
WIND  AT  FLIGHT  ALTITUDE. 
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FIG.  2.  WINDS  REQUIRED  FOR  'CUT-OFF1  OF  SONIC  BOOM. 


STANDARD  ATMOSPHERE 
ALTITUDE  =  50,000  ft. 

A,  M  =  1.3,  W  =  64.5  knots 

B,  M  =  1.6,  W  =  64.5  knots 

C,  M  =  2.05,  W  =  64.5  knots 

D,  M  =  1.5,  W  =  200  knots 

E,  M  =  1.6,  W  =  200  knots 
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1.  No  wind 

2.  Tail  wind 

3.  Side  wind 

4.  Head  wind 


Flight 


FIG.  3.  EFFECTS  OF  MACH  NUMBER  AND  WIND  MAGNITUDE  ON 
SONIC  BOOM  CORRIDOR. 
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FIG.  6.  COMPARISON  OF  PRESENT  RESULTS  WITH  PUBLISHED  DATA. 


Corridor  Widths,  mis 


160 

120 


Altitude 

M 

W 


50,000  ft. 
2.05 

64.5  knots 


80 


CORRIDORS  REPRESENT: 


40 

0 


1. 

2. 

3. 

4. 


No  wind 
Tail  wind 
Side  wind 
Head  wind 


Flight 
Dir 'n. 


e 

o 

£ 

o 

-H 

c 

o 

to 


40t- 


Std.  Atm 

AT  =  Tg-Th  =  128.70 


°F 


80 


Non-Std.  Atm.  (Winter) 
AT  =  Tg-Th  =  77.76  °F 


120 


1601- 


fig.  7.  EFFECTS  OF  GROUND  TEMPERATURE  ON  SONIC  BOOM  CORRIDOR 


